Infrared probe of the anomalous magnetotransport of highly oriented pyrolytic 

graphite in the extreme quantum limit 



O 

o 

> 
O 

o 



X3 

o 



> 

(N 



o 



I 

o 
o 



X 
S3 



Z.Q. Li,i'B S.-W. Tsai,2 W.J. Padilla,i'[l] S.V. Dordevic,^ K.S. Burch,i'0 Y.J. Wang,"* and D.N. Basoyi 

^Department of Physics, University of California, San Diego, La Jolla, California 92093, USA 
^Department of Physics, University of California, Riverside, California 92521, USA 
'^Department of Physics, The University of Akron, Akron, Ohio 44-325, USA 
'''National High Magnetic Field Laboratory, Tallahassee, Florida 32310, USA 
(Dated: February 4, 2008) 

We present a systematic investigation of the magneto-reflectance R{uJ, B) of highly oriented py- 
rolytic graphite (HOPG) in magnetic field B up to 18T. From these measurements, we report the 
determination of lifetimes r associated with the lowest Landau levels in the quantum limit. We find 
a linear field dependence for inverse lifetime 1/t{B) of the lowest Landau levels, which is consis- 
tent with the hypothesis of a three-dimensional (3D) to ID crossover in an anisotropic 3D metal 
in the quantum limit. This enigmatic result uncovers the origin of the anomalous linear in-plane 
magneto-resistance observed both in bulk graphite and recently in mesoscopic graphite samples. 



L INTRODUCTION 

Intense magnetic fields can radically alter properties 
of materials and produce a myriad of novel effects and 
novel states of matter.-i^ One intriguing example is the 
field-induced transformation expected in an anisotropic 
three-dimensional (3D) metal. Provided only the lowest 
Landau levels (LLs) are populated (a condition referred 
to as the quantum limit), coherent charge motion in the 
direction perpendicular to the magnetic field is arrested. 
Therefore magnetic field initiates a crossover from 3D to 
ID transport along the direction of the field4 A suitable 
candidate system to exhibit this phenomenon is graphite. 
In this anisotropic semimetal comprised of weakly cou- 
pled graphene sheets, a combination of low carrier den- 
sity and small effective masses allows one to achieve the 
quantum limit condition in fairly modest magnetic fields 
along the c-axis (T-ST).*^ A manifestation of the 3D to 
ID crossover is the field dependence of the relaxation 
rate 1/r associated with the lowest LLsj^ This funda- 
mental electronic characteristic is difficult to infer from 
transport studies but can be directly probed via infrared 
(IR) spectroscopy. 

Here, via IR spectroscopy of graphite we probe life- 
times T associated with the lowest LLs in the quantum 
limit with field perpendicular to the graphene planes. We 
find that 1/t oc B for the lowest LLs, which is distinct 
from power laws typically found in other classes of low 
dimensional systems^ but consistent with the hypoth- 
esis of a 3D to ID crossover^ Empowered by the ex- 
perimental data for lifetimes, we were able to elucidate 
the anomalous linear magneto-resistance in samples of 
weakly interacting graphene sheets. This peculiar form 
of magneto-resistance was discovered nearly 40 years ago 
in bulk graphite^ ^ and recently observed in mesoscopic 
samples comprised of only few graphene layers^^i 

II. EXPERIMENTAL METHODS 

HOPG samples studied in this work were ex- 
tensively characterized through magneto-transport 



measurementsi^ They have ab-plane surfaces with 
typical dimensions 6x6 mm^. The zero field reflectance 
R{iLj) was measured over broad frequency (15 to 25,000 
cm^^) and temperature (10 to 292 K) ranges at UCSD. 
The complex optical conductivity a{uj) — ai{uj) + ia2{io) 
in zero field was obtained from Kramers-Kronig (KK) 
transformation of the refiectance data. The Hagen- 
Rubens formula was employed in the low lo region. 
KK-derived results are consistent with those obtained 
by ellipsometry. 

The magneto-refiectance R{w, B) were measured in the 
Faraday geometry: E vector in the graphene plane and 
magnetic field B parallel to the c axis of the sample. 
Our in-house apparatus enables absolute measurements 
oiR{uj,B) in fields up to 9T over a broad frequency range 
(15 to 4,000 cm^^).^ Data in selected magnetic fields were 
acquired using a sub-THz Martin-Puplett interferometer 
and extend down to 7 cm^^. The magneto- reflectance 
ratios R{uj , B) / Rluj , B = 0) in the range 20-3,000 cm"! 
with B=7-18 T were obtained at the National High Mag- 
netic Field Laboratory in Tallahassee. All data in mag- 
netic fleld were taken at 5K. 

HI. ZERO FIELD DATA 

The zero fleld reflectance -R(w) and dissipative part 
of the optical conductivity ai{Lo) are plotted in Fig. 1. 
The room temperature i?(w) spectrum is consistent with 
earlier studies in the overlapping frequency rangeJ^ The 
R{uj) spectrum shows a typical metallic behavior: a grad- 
ual depression of reflectance with increasing frequency 
towards a plasma minimum at about 24,000 cm~^, fol- 
lowed by a sharp onset due to the interband transition 
between tt bands at the M points of the Brillouin zonei^ 
An edge-like structure in the far-IR reflectance shows an 
interesting temperature dependence with softening of the 
edge from 400 cm^^ at 292 K down to 200 cm^^ at 10 K. 
The zero fleld (Ti{lo) spectra are dominated by the Drude 
component but in addition show a substantial contribu- 
tion throughout the mid-IR range. This is detailed in 
the bottom inset of Fig. 1 displaying a broad resonance 
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corresponding to the edge observed in raw R(uj) data. As 
temperature is reduced the Drude mode narrowed indica- 
tive of strong depression of the quasiparticle scattering 
rate 1/t. We also found that mid-IR resonance softens 
with the center peak shifting from 1,000 cm~^ at 292 K 
down to 400 cm^^at 10 K. 

We assign the mid-IR resonance to the interband tran- 
sitions between the tt bands near K and H points of 
the Brillouin zone,— as shown in Fig. 2. The conven- 
tional Slonzcewski-Weiss-McClure (SWMC) band theory 
of graphite^ prescribes the following relations for the 
onsets of interband transitions near K and H pointsjii 
AEk < 2{Ep - 272) and I^Eh = A - 2Ef, where A and 
72 are the band parameters. Therefore, the lowering of 
the Fermi energjs^i^ and tt bands near the K poini*^ with 
the increase of T is expected to cause hardening of the 
transitions at K and H points, as shown in Fig. 2. The 
observed trends are consistent with the depression of Ep 
by 40% at 292K compared to the value at 10 K assuming 
accepted values for the Fermi energy and the band pa- 
rameters (at lOK): Ef = -25meV, A = 0.04 ~ QMeV^ 
and 72 = — 0.02eT^. Although this effect has not been 
directly observed before, the evolution of the resonance 
structure in ai{oj) with T is consistent with the well es- 
tablished notion that the dominant contribution to the 
carrier dynamics and electrodynamics of graphite is asso- 
ciated with K and H points of the Brillouin zonc.^ There- 
fore, our zero field results are consistent with the SWMC 
band model. The applicability of SWMC model was cor- 
roborated by recent angle resolved photoemission spec- 
troscopy studies of HOPC— 

IV. MAGNETO-REFLECTANCE DATA 

Insights into the anomalous magneto-transport can 
be provided by a systematic investigation of magneto- 
reflectance R{oj, B) spectra, which are displayed in a w/B 
scaled plot in Fig. 3. To highlight field-induced changes 
of the optical properties, at each field B the zero field re- 
flectance spectrum B = 0) is also scaled by the same 
factor B and shown together with B) for compar- 
ison. Two new features are observed in the high-field 
spectra: i) the reflectance at low frequency is strongly 
suppressed and no longer extrapolates to unity in the 
limit of ^ 0; ii) the R(uj,B) spectra reveal a series of 
field-dependent resonances extending up to 3,000 cm~^ 
(see also the bottom panel of Fig. 4). The dramatic sup- 
pression of the low frequency reflectance R{u} — > 0, B) is 
an IR counterpart of the enormous positive magnetore- 
sistance of graphite exceeding 10'^ in 18 T field. ^ The 
resonance features in i?(aj, B) originate from the inter- 
LL transitions. An inspection of Fig. 3 shows that the 
frequencies of the strongest resonances reveal a linear de- 
pendence on magnetic field. 

The rather complicated form of the i?(w, B) spectra 
and its systematic evolution with magnetic field are both 
in accord with the SWMC band theory of HOPG. To 
show this we analyzed R{uj, B) spectra displayed in Fig. 3 



using a magneto-optical Lorentzian model for the com- 
plex conductivity due to inter-LL transitions a^{ijj)'i^ 

= Vi^^_ (1) 

This equation yields a series of resonance peaks at LOj+ 
and bjj- at the poles of the left- and right-hand circu- 
larly polarized conductivity (w) . Each peak is defined 
through the linewidth 1/tj± and the oscillator strength 
'^'^j± ■ Inter-LL transitions at both K and H points of the 
Brillouin zonei^ii are included in our fitting. According 
to SWMC theory, more than a dozen inter-LL transi- 
tions are likely to contribute to i?((j-', B) in the studied 
frequency range. However, the parameters in Eq. (1) are 
strictly constrained by both the analytical properties of 
(j±{u})'^^ and the selection rules of inter-LL transitions.^ — 
Following these constraints, Eq. (1) is used to construct 
a+{ijj) and <J-{lo)^ which are then used to evaluate the 
reflectance spectra R{uj) = \[R+{lo) + i?_(a;)]. In addi- 
tion, several fleld-independent Lorentzian oscillators are 
employed in both (Ji+{ijj) and cri_(a;) to account for the 
absorption above 3,000 cm~^, as displayed by the green 
(light gray) curve in Fig. 4. We observed an additional 
resonance feature in R{lo^ B) below 8.5T in the range 50- 
100 cm^^/T in Fig. 3, which is an expected consequence 
of the quantum limit. The resonances in R(u}, B) below 
and above 8.5T are roughly at the same positions in the 
Lj / B plot, with small deviations in the high frequency 
resonances. The two different sets of magneto-optical os- 
cillators at uj.j±/ B employed to flt experimental R{lo, B) 
are summarized in table 1. We use the transition ener- 
gies from refiii to account for the inter-LL transitions at 
the H point. These transitions are at much higher en- 
ergy compared to those at the K point and constitute a 
broad background in R{uj,B) spectra.-i^ We stress that 
all fleld-dependent oscillators employed in our model can 
be assigned to the inter-LL transitions at K and H points 
within the conventional SWMC band model. ^ 

The model spectra at representative fields obtained 
from the above analysis are depicted in Fig. 3 and 
4. Both the peak structures and their field dependence 
in experimental R(uj,B) spectra are reproduced by our 
analysis. Of special interest is the lowest inter-LL transi- 
tion (LL-l)-(LLO) shown in Fig. 5. This transition is well 
separated from higher energy resonances, therefore the 
reflectance in the vicinity of this mode is hardly affected 
by other higher energy modes in our model, which allows 
us to extract the parameters for this particular transi- 
tion from the above analysis. Moreover, in the quantum 
limit LL-1 and LLO are the only occupied levels;^ which 
are directly related to the magneto-transport properties. 
Therefore, special attention will be given to the (LL-1)- 
(LLO) transition. Both the resonance frequency of this 
transition and its linewidth 1/t follow a linear field de- 
pendence as shown in Fig. 5. This is evident from a close 
inspection of the experimental R{u!, B) spectra in Fig. 3. 
The power law of 1/t{B) below 6T field is yet to be ex- 
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plored, because at low fields the (LL-l)-(LLO) transition 
occurs at low energies beyond our detection limit. 



V. DISCUSSION 

Our magneto-reflectance data are in accord with the 
conventional SWMC band model of graphitci^ In the 
SWMC model, the inter-LL transitions at K and H points 
have different field dependence: transitions at the K 
point reveal "normal" quasiparticle behavior (linear B 
scaling);^ whereas transitions at the H point show Dirac 
behavior at high fields (~ Vb scaling). -'^^ Inter-LL transi- 
tions at both K and H points of the the Brillouin zone are 
needed to fully account for our results on HOPG, which 
is consistent with previous magneto-optical study of bulk 
graphite samples.^'' Recently, inter-LL transitions with 
^/B scaling were observed in magneto-transmission study 
of ultrathin graphite layerSfi^ which may come from tran- 
sitions at the H point due to inter-layer coupling between 
graphene sheets. 

We now show that the observed 1/t on B dependence 
for the (LL-l)-(LLO) transition (Fig. 5) is consistent 
with the notion of 3D to ID crossover in an anisotropic 
3D metal in the quantum limit4 The 1/t (x B behav- 
ior in Fig. 5 most likely originates from the linear field 
dependence of the broadening of both LL-1 and LLO,^° 
which defines the scattering rate of the quasiparticlcs as- 
sociated with these levels. In the quantum limit, the 
carriers have finite momentum only along the field di- 
rection (c-axis) and move diffusively in the transverse 
direction (ab plane) when scattered by impurities. The 
carrier-impurity scattering rate 1 /r in the quantum limit 
is given by 1/t{B) = 2T:nQUQVBr where the density of 
neutral impurities ng and impurity potential uq are B- 
independent. The 3D density of states vb increases lin- 
early with B in the quantum limit in graphite,^ giving 
rise to the linear B dependence of 1/t. The increase of 
1/r by a factor of 3 from 6T to 18T is consistent with 
a threefold increase of vb in this range.— The interac- 
tions between carriers can give additional corrections to 
1/t{B), which have weak logarithmic field dependence. 
Therefore our data for the field dependence of the relax- 
ation rate are in accord with the specific prediction of a 
model involving 3D to ID crossover in the quantum limit. 

Based on the observed 1/t oc B result, we propose an 
explanation for the anomalous linear magneto-resistance 
found both in mesoscopic graphite systems! and in bulk 
graphite?^ First, the giant positive magneto-resistance 
of graphite is clearly related to the electronic spectral 
weight transfer of the Drude mode in the B=0 spectrum 
to finite energies in magnetic field (Fig. 4). The spec- 
tral weight transfer results in a drastic suppression of the 
conductivity in the uj —^ limit by more than 10^. Data 
in Figs. 3, 4 furthermore show that the lowest energy 
electronic excitation in the quantum limit is the (LL-1)- 
(LLO) transition revealing the linear field dependence of 
the relaxation rate. Let us explore the consequences of 



these findings for magneto-resistance that we will analyze 
following the standard form appropriate for the quantum 
limit as shown ini22d 

[p^ + P2Y + {Ri + R2?B^ ^' 



where p^.^. is the ab-plane resistivity, p^ and =l/niqj 
are the resistivity and Hall coefficient of the electrons 
(i=l) and holes (i=2), and n^ and q^ are the density and 
charge of the carriers. Note that in the quantum limit 
the off-diagonal term in the resistivity tensor can still be 
expressed as RiB as shown in Ref4. We emphasize that 
here Eq. (2) follows from a fully quantum mechanical 
treatment yielding the magnetic field dependence of p.^ 
and Ri. This is distinct from the semiclassical model of 
magnetoreistance, in which p^ and Ri are constants. At 
high fields, the B^ term dominates over the zero field 
term in the numerator, therefore the latter can be ne- 
glected. In the quantum limit, transport in the ab plane 
is diffusive with p^ given by ~ e^usl^r^ where D 
~ ^b/t is the diffusion coefficient and Ib = 1/VeB is 
the magnetic length. From vb oc B^ and D ^ I'b/t 
oc B/B ^1 {1/t cx B from our results), we obtain 
1/ Pi oc B. Because the net carrier density in the sample 
7ii — is constant at high field^^'^ and ni « n2 oc _B, 
we obtain {Ri + i?2)B = ^"'^^""^ ^ oc 1/B. This expres- 
sion is valid in the magnetic field regime where graphite 
shows linear magneto-resistance. In even higher field, 
7T-1 — 712 decreases with field due to the so-called mag- 
netic freeze-out effect, and correspondingly a satura- 
tion of magneto-resistance is observed.^^ Combined with 
Ri = l/rije oc 1/B, we find p^^ oc B. Therefore, 
our study suggests that the anomalous linear magneto- 
resistance of graphite originates from quantum transport 
in magnetic field. The above analysis can also be used to 
account for the linear magneto-resistance in mesoscopic 
graphite-based electronic systems^ where interlayer cou- 
pling between graphene sheets is important. 

The use of a two-band model in our analysis is corrob- 
orated by a recent study of mesoscopic graphite field- 
effect transistors)! which shows that the slope of the 
magneto-resistance can be tuned by changing the car- 
rier compensation n,, — Uh consistent with Eq (2). Pre- 
viously, a model based on magnetic-field-dependent scat- 
tering range of ionized impurity scattering^ was proposed 
to explain the linear magneto-resistance in graphite. Due 
to the field-dependent screening of charged impurities, 
the scattering rate in this model has a field dependence 
of 1/t ~ B/{c + B)'^, where c is a constant.^ This predic- 
tion is in contrast to our observation 1/r oc B. There- 
fore, our study shows that impurities in graphite are 
likely to be uncharged and the ionized impurity scat- 
tering model^ may not be applicable. Indeed, previous 
magneto-transport studies have shown the limitations of 
the ionized impurity model i^^^ Recently, an interesting 
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model based on Dirac fermions^^ is proposed to account 
for linear magneto-resistance in graphite. However, our 
studies of HOPG have clearly demonstrated that Dirac 
fermions do not play a dominant role in the magneto- 
transport in graphite. 

VI. SUMMARY 

In this study, we carried out systematic IR magneto- 
optical experiments on HOPG. We show that the opti- 
cal properties of HOPG in both zero field and magnetic 
fields can be understood within the conventional SWMC 
band theory of graphite.^ We observed an unconventional 
linear field dependence of the inverse lifetime l/r asso- 
ciated with the lowest LLs. This observation is consis- 
tent with the hypothesis of 3D to ID crossover in an 



anisotropic 3D metal in the quantum limit. ^ This new re- 
sult has allowed us to address the origin of the anomalous 
linear magneto-resistance in mesoscopic graphite-based 
electronic systems^ as well as in bulk graphite.^ Our work 
has demonstrated the potential of IR spectroscopy for the 
study of LL dynamics and novel magneto-transport phe- 
nomena in carbon based materials such as mesoscopic 
graphite and graphene^^Z^ which are directly related 
to the novel functionalities of magneto-electronic devices 
based on these materials. 
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TABLE I: The magneto-optical Lorentzian oscillators used in the analysis of Ti{uj, B) and their possible assignments. uj+I B and 
u^/B (in cm~^/r) are frequencies (normalized by magnetic field B) of the oscillators in o-+(tj) and a_(cij), repectively, used to 
account for the inter-Landau-level (LL) transitions at K-point. Two sets of oscillators are used in two field ranges: (h) B >8.5T 
and {I) B ^8.5T. The theoretical transition frequencies are labeled following Nakao's notationi^: hn-em refers to the transition 
from the n-th hole LL to the m-th electron LL, etc. The transitions marked with stars are allowed only below the quantum 
limit. The last 8 columns list u>pj± and 1/Tj± (both in cm^^) used for the model spectrum R(t<;) at 16 and 8.5T. The parameters 
of the oscillator used for m^=0,-l transitions at H-point are: 16 T: loq — 1212 cm~^,ajp = 4243cm-S 1/r = 500cm-^ 8.5T: 
Wo = 818 cm"\a;p = 3211cm"\l/r = 300cm"\ 
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FIG. 1: (Color online) Zero field ab-plane reflectance R{lj) 
(top panel) and real part of the optical conductivity g\{ijj) 
(bottom panel) of graphite. The thin dashed lines in the 
bottom panel are Drude fits of the ai{uj) spectra in the low 
frequency region. The bottom inset displays the absorption 
peak of (Ti (u) in the mid-IR region, which shifts to high energy 
with increasing temperature. 



FIG. 2: (Color online) Schematics of the band structure of 
graphite at K and H points of the Brillouin zone at lOK and 
292K. The arrows indicate the interband transtions observed 
in the bottom inset of Fig. 1. 



FIG. 3: (Color online) The uj/B scaled plot of the ab-plane 
reflectance spectra R(uj, B) in several magnetic fields. At each 
field B, the zero field refiectance spectrum R{uj, i? = 0) is also 
scaled by the corresponding factor B for comparison. Every 
set of spectra is offset by 30% for clarity. The model spectra 
are obtained from the analysis detailed in the text. 
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FIG. 4: (Color online) Top panel: (Ti(tij) at zero field and 
model spectra (Ji+(tj) and ai-{uj) at 16T obtained from the 
analysis detailed in the text. The black square on the left axis 
represents the DC conductivity aoc at OT. The green (light 
gray) curve is the high frequency field-independent absorp- 
tion in both ai-{(jj) and cti+(lj). Bottom panel: experimen- 
tal R(tj) spectra at OT and 16T and model spectra R(c<j) at 
16T. The arrows indicate the resonance features in R(ci;) and 
ai-{uj) due to the (LL-l)-(LLO) transition. 



FIG. 5: (Color online) Resonance frequency (left panel) and 
linewidth (middle panel) of the (LL-l)-(LLO) transition in 
magnetic field. The dashed curves in these two panels are 
linear fits to the data. The dotted curves in the left panel are 
u! ~ y/B fits. Right panel: a schematic of the LL-1 and LLO 
Landau levels in the quantum limit (adapted from Ref.-^). 
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